3.1 State Space Models 

In this section we study state space models of continuous-time lin¬ 
ear systems. The corresponding results for discrete-time systems, 
obtained via duality with the continuous-time models, are given in 
Section 3.3. 

The state space model of a continuous-time dynamic system 
can be derived either from the system model given in the time 
domain by a differential equation or from its transfer function 
representation. Both cases will be considered in this section. Four 
state space forms—the phase variable form (controller form), the 
observer form, the modal form, and the Jordan form—which are 
often used in modern control theory and practice, are presented. 


3.1.1 The State Space Model and Differential Equations 

Consider a general nth-order model of a dynamic system repre¬ 
sented by an nth-order differential equation 


<Fy{t) , d^-^y{t) , 

+ CLn-l -TTI-n-h 


dt 


n 


dt^-^ 


+ + my{t) 


dt 


, d^uit) , cP ^u(t) 

= + K-l ^ + 


dt 


dt 


+ + bou(t) 

at 


(3.1) 

At this point we assume that all initial conditions for the above 
differential equation, i.e. y(()~),dy(0~)/dt, ...,(r'~^y(()~)/dt^~^, 
are equal to zero. We will show later how to take into account the 
effect of initial conditions. 


In order to derive a systematic procedure that transforms a 
differential equation of order n to a state space form representing 
a system of n first-order differential equations, we first start with 
a simplified version of (3.5), namely we study the case when no 
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derivatives with respect to the input are present 

(Pyit) (r~^y{t) dyit) 

—;—- + an-i —;-+ • • • + ai —+ CLoy(t) = u(t) (3.2) 

Introduce the following (easy to remember) change of variables 

xi{t) = y{t) 

«») - ^ 


xsit) 


(fy(t) 


(3.3) 


which after taking derivatives leads to 


= r.(t) 

dt ^ dt 
dx 2 {t) . cfy(t) , 


dx2 (t) 
dt 

dxs (t) 
dt 


= X2 = 


X3(t) 


d^y(t) .. 

= X3 = ^^3 = X4{t) 


dxnit) _ . _ d!^y{t) 

dt ^ dt^ 

. . dy(t) dl^yit) d^~^y(t) 

’ dt dt^ dt^-^ ^ ’ 

= —aoxi{t) — aiX2{t) — ... — a2X3{t) — ... — an-iXn{t) + u{t) 


(3.4) 
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The state space form of (3.8) is given by 



(3.5) 


with the corresponding output equation obtained from (3.7) as 

xi(t) 

X2(t) 

y(t) = [l 0 ••• 0] ! (3.6) 

Xn-l(t) 

_ Xnit) 

The state space form (3.9) and (3.10) is known in the literature as 
the phase variable canonical form. 

In order to extend this technique to the general case defined 
by (3.5), which includes derivatives with respect to the input, we 
form an auxiliary differential equation of (3.5) having the form of 
(3.6) as 

+ ^n-l yyn-l 1“ ' ' ' 3“ ^ = u(t) (3.7) 


dt^ 
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for which the change of variables (3.7) is applicable 

xi{t) = ^(t) 
d^(t) 


X2(t) = 
xsit) 


dt 
dt^ 


(3.8) 


Xn(t) 


d^-^m 

dt^-^ 


and then apply the superposition principle to (3.5) and (3.11). Since 
^(t) is the response of (3.11), then by the superposition property 
the response of (3.5) is given by 


y(t) = bom + hi^ + + 


dt 


dt^ 


+ h 


d'm 


n 


dt 


n 


(3.9) 


Equations (3.12) produce the state space equations in the form 
already given by (3.9). The output equation can be obtained by 
eliminating (F^it)/dt^ from (3.13), by using (3.11), that is 


d^m 

dt^ 


u(t) — an-iXn(t) — ... — aiX 2 (t) — aoxi(t) 


This leads to the output equation 


y(t) = [(^0 - CLobn) (bi - aibn) 



Xl(t) 

X 2 (t) 



+bnu(t) 

(3.10) 

It is interesting to point out that for bn = 0, which is almost always 
the case, the output equation also has an easy-to-remember form 
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given by 

xi(t) 

y(t) = [bo bi ••• bn-i] (3.11) 

_Xn(t)_ 

Thus, in summary, for a given dynamic system modeled by dif¬ 
ferential equation (3.5), one is able to write immediately its state 
space form, given by (3.9) and (3.15), just by identifying coeffi¬ 
cients at and 6*, i = 0,1,2,n — 1, and using them to form the 
corresponding entries in matrices A and C. 

Example 3.1: Consider a dynamical system represented by the 
following differential equation 

+ 6^^^- — 2^^^- + — 5^^^- + 3^ = + Au 

where stands for the zth derivative, i.e. y^^'^ = dy/dt\ 
According to (3.9) and (3.14), the state space model of the above 
system is described by the following matrices 


" 0 1 0 0 0 0 1 To" 

0 0 1 0 0 0 0 

000100 T3_0 

000010 ’ 0 
0 0 0 0 0 1 0 

-3 5 -1 0 2 -6] [l 


C = [4 1 0 7 0 0], D = 0 


o 
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3.1.2 State Space Variables from Transfer Functions 

In this section, we present two methods, known as direct and 
parallel programming techniques, which can be used for obtaining 
state space models from system transfer functions. For simplicity, 
like in the previous subsection, we consider only single-input 
single-output systems. 

The resulting state space models may or may not contain all the 
modes of the original transfer function, where by transfer function 
modes we mean poles of the original transfer function (before 
zero-pole cancellation, if any, takes place). If some zeros and 
poles in the transfer function are cancelled, then the resulting state 
space model will be of reduced order and the corresponding modes 
will not appear in the state space model. This problem of system 
reducibility will be addressed in detail in Chapter 5 after we have 
introduced the system controllability and observability concepts. 

In the following, we first use direct programming techniques to 
derive the state space forms known as the controller canonical form 
and the observer canonical form; then, by the method of parallel 
programing, the state space forms known as modal canonical form 
and Jordan canonical form are obtained. 

The Direct Programming Technique and Controller Canon¬ 
ical Form 

This technique is convenient in the case when the plant transfer 
function is given in a nonfactorized polynomial form 

v(5) _ bnS^ + bn-lS^ ^ + * * * + bis + 6o ^ 

U (5) + aji—is^ ^ 5 + uq 

For this system an auxiliary variable V{s) is introduced such that 
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V{s) 

U{s) 

y{s) 

V{s) 



101 

is split as 


1 

(3.13a) 

+ a-n—lS^ ^ 5 + Uq 

i>nS^ + bn—lS^ ^ 5 + 6o 

(3.13b) 


The block diagram for this decomposition is given in Figure 3.1. 



Figure 3.1: Block diagram representation for (3.17) 

Equation (3.17a) has the same structure as (3.6), after the 
Laplace transformation is applied, which directly produces the state 
space system equation identical to (3.9). It remains to find matrices 
for the output equation (3.2). Equation (3.17b) can be rewritten as 

y(s) = b„s"y(s)+b„_is’‘-V(s)+- • •+6isy(s)+6oy(s) (3.14) 

indicating that y(t) is just a superposition of v(t) and its derivatives. 
Note that (3.17) may be considered as a differential equation in the 
operator form for zero initial conditions, where s = d/dt. In that 
case, V(s), Y(s), and U(s) are simply replaced with v(t), y{t), 
and u(t), respectively. 

The common procedure for obtaining state space models from 
transfer functions is performed with help of the so-called transfer 
function simulation diagrams. In the case of continuous-time 
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systems, the simulation diagrams are elementary analog computers 
that solve differential equations describing systems dynamics. They 
are composed of integrators, adders, subtracters, and multipliers, 
which are physically realized by using operational amplifiers. In 
addition, function generators are used to generate input signals. The 
number of integrators in a simulation diagram is equal to the order 
of the differential equation under consideration. It is relatively easy 
to draw (design) a simulation diagram. There are many ways to 
draw a simulation diagram for a given dynamic system, and there 
are also many ways to obtain the state space form from the given 
simulation diagram. 

The simulation diagram for the system (3.17) can be ob¬ 
tained by direct programming technique as follows. Take n 
integrators in cascade and denote their inputs, respectively, by 

Use formula (3.18) to construct 
y(t), i.e. multiply the corresponding inputs to integrators 

by the corresponding coefficients b^ and add them using an adder 
(see the top half of Figure 3.2, where 1/s represents the integrator 
block). From (3.17a) we have that 


= u{t) — an-iv^^ — • • • — aiv^^\t) — aov(t) 

which can be physically realized by using the corresponding feed¬ 
back loops in the simulation diagram and adding them as shown 
in the bottom half of Figure 3.2. 
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Figure 3.2: Simulation diagram for the direct 
programming technique (controller canonical form) 

A systematic procedure to obtain the state space form from 
a simulation diagram is to choose the outputs of integrators as 
state variables. Using this convention, the state space model for 
the simulation diagram presented in Figure 3.2 is obtained in a 
straightforward way by reading and recording information from 
the simulation diagram, which leads to 



"o 1 0 . o' 


1 

o 

1_ 

0 0 1 0 ••• 0 


0 

o ••• 

x(0 + 

• 

0 0 0 ••• 0 1 


0 

—no —oi —02 . —On-l 


1 

_1 



(3.15) 
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and 


y(f) = K^O - (^1 - aihn) 

+hnu{t) 


(^n—1 ]x(t) 

(3.16) 


This form of the system model is called the controller canon¬ 
ical form. It is identical to the one obtained in the previous sec¬ 
tion—equations (3.9) and (3.14). Controller canonical form plays 
an important role in control theory since it represents the so-called 
controllable system. System controllability is one of the main con¬ 
cepts of modern control theory. It will be studied in detail in 
Chapter 5. 

It is important to point out that there are many state space forms 
for a given dynamical system, and that all of them are related by 
linear transformations. More about this fact, together with the 
development of other important state space canonical forms, can 
be found in Kailath (1980; see also similarity transformation in 
Section 3.4). 

Note that the MATLAB function tf2ss produces the state 
space form for a given transfer function, in fact, it produces the 
controller canonical form. 


Example 3.2: The transfer function of the flexible beam from 
Section 2.6 is given by 


, . 1.65s^ - O-SSls"® - 576 s 2 + 90.6s + 19080 

^ + 0.996s5 + 463 s 4 + 97.8s3 + 12131s2 + 8.11s 


Using the direct programming technique with formulas (3.19) and 
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(3.20), the state space controller canonical form is given by 


X = 


and 


"0 

1 

0 

0 

0 

0 


"0" 

0 

0 

1 

0 

0 

0 


0 

0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

0 

1 

0 

X+ 

0 

0 

0 

0 

0 

0 

1 


0 

0 

-8.11 - 

12131 

-97.8 

-463 

-0.996 


1 

y 

= [19080 

90.6 

-576 

-0.331 

1.65 0]x 



o 

Direct Programming Technique and Observer Canonical 
Form 

In addition to controller canonical form, observer canonical 
form is related to another important concept of modern control 
theory: system observability. Observer canonical form has a very 
simple structure and represents an observable system. The concept 
of linear system observability will be considered thoroughly in 
Chapter 5. 

Observer canonical form can be derived as follows. Equation 
(3.16) is written in the form 

Y(s) ( 5 ^ + On-is^ ^ + • • • + + no) 

= U(s)(b„s’‘+b„-is’‘-^ + ---+bis + bo) '■ ' 

and expressed as 

1^(5) = — -(^On-is^ ^ + • • • + + no)^(' 5 ) 

s 

+ —( 5 ) (bnS^ + bn-lS^ ^ + * ’ ’ + ^> 1-5 + b^j 

^ (3.18) 
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leading to 


y(s) = -ia„-iy(s) - \an- 2 Y(s) - 


1 


aiy(s) 


-4aoy(s) + hnU{s) + h„-lU{s) + \bn- 2 U(s)+ 

^ 11 o 


■■■ + J=jhU(s) + j^l^U{s) 

(3.19) 

This relationship can be implemented by using a simulation di¬ 
agram composed of n integrators in a cascade, and letting the 
corresponding signals to pass through the specified number of in¬ 
tegrators. For example, terms containing 1/s should pass through 
only one integrator, signals an- 2 y{t) and bn- 2 u{t) should pass 
through two integrators, and so on. Finally, signals aQy(t) and 
bou(t) should be integrated n-times, i.e. they must pass through 
all n integrators. The corresponding simulation diagram is given 
in Figure 3.3. 
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Figure 3.3: Simulation diagram for observer canonical form 


Defining the state variables as the outputs of integrators, and 
recording relationships among state variables and the system out¬ 
put, we get from the above figure 

y{t) = Xn{t) + hnU{t) 

Xi{t) = -a{)y{t) + hu{t) = -a{)Xn{t) + (6o - aobn)u(t) 

X2(t) = -aiy(t) + biu(t) + xi(t) 

= xi(t) - aiXn(t) + (bi - aibn)u(t) 

xsit) = -a2y(t) + b2u(t) + X2(t) 

= X2(t) - a2Xn(t) + (62 - a2bn)u(t) 

Xn(t) = -an-iy(t) + bn-iu(t) + Xn-l(t) 

= + (bn-1 - an-lbn)u(t) 

The matrix form of observer canonical form is easily obtained 
from (3.24) and (3.25) as 


(3.20) 


(3.21) 
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and 


"o 

0 ... 

... 0 

-no 


bo- 

- aobn 

1 

0 ... 

... 0 

—ai 


bi- 

- aibn 

0 

1 

0 1 

... : 

— <32 

x(0+ 

b2- 

- CL2bn 

• 

• • , 

0 

<3n—2 



• 

0 

0 ... 

0 1 

<3n—1 


bn—l 

“ O^n—lbn 


(3.22) 

^(t) = [0 ••• 0 1 ]x(t) + (3.23) 


Example 3.3: The observer canonical form for the flexible 
beam from Example 3.2 is given by 


"o 

0 

0 

0 

0 

0 


" 19080 " 

1 

0 

0 

0 

0 

-8.11 


90.6 

0 

1 

0 

0 

0 

-12131 

x+ 

-576 

0 

0 

1 

0 

0 

-97.8 

-0.331 

0 

0 

0 

1 

0 

-463 


1.65 

0 

0 

0 

0 

1 

-0.996 


0 


?/ = [0 0 0 0 0 l]x 


Observer canonical form is very useful for computer simulation 
of linear dynamical systems since it allows the effect of the sys¬ 
tem initial conditions to be taken into account. Thus, this form 
represents an observable system, in the sense to be deflned in 
Chapter 5, which means that all state variables have an impact 
on the system output, and vice versa, that from the system output 
and state equations one is able to reconstruct the state variables 
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at any time instant, and of course at zero, and thus, determine 

in terms of the original initial conditions 
y(0~),dy(0~)/dt, ...,(r~^y(0~)/dt^~^. For more details see Sec¬ 
tion 5.5 for a subtopic on the observability role in analog computer 
simulation. 

Parallel Programming Technique 

For this technique we distinguish two cases: distinct real roots 
and multiple real roots of the system transfer function denominator. 

Distinct Real Roots 

This state space form is very convenient for applications. 
Derivation of this type of the model starts with the transfer func¬ 
tion in the partial fraction expansion form. Let us assume, without 
loss of generality, that the polynomial in the numerator has degree 
of m < n, then 


U{s) 


_ Pm{s) _ 

(S + pi)(s + P2) "'{S+Pn) 

kl ^ k 2 ^ ^ kn 

S-\-Pi S-\-p 2 S+Pn 


(3.24) 


Here pi,p 2 ,---,Pn are distinct real roots (poles) of the transfer 
function denominator. 

The simulation diagram of such a form is shown in Figure 3.4. 
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Figure 3.4: The simulation diagram for the parallel 
programming technique (modal canonical form) 

The state space model derived from this simulation diagram is 
given by 

—pi 0 . 0 1 

0 —P2 ■•.••• 0 1 

x(t) = : ■ •. ■ •. ■ •. : x(t) + : u{t) 

\ 0 : (3.25) 

0 0 • • • 0 —Pn 1 

y(t) = [ki k 2 . kn]^{t) 

This form is known in the literature as the modal canonical form 
(also known as uncoupled form). 

Example 3.4: Find the state space model of a system described 
by the transfer function 

y(5) _ (5 + 5)(5 + 4) 

U ( 5 ) (5 + 1)(5 + 2)(5 + 3 ) 
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using both direct and parallel programming techniques. 

The nonfactorized transfer function is 

Y (s) 5^ + 95 + 20 

U(s) 5^ + 65^ + 115 + 6 

and the state space form obtained by using (3.19) and (3.20) of the 
direct programming technique is 


" 0 

1 

0 ■ 


"o' 

0 

0 

1 

x + 

0 

-6 

-11 

-6 


1 


?/ = [20 9 l]x 

Note that the MATLAB function tf2ss produces 


'-6 

-11 

-6" 


"l" 

1 

0 

0 

x + 

0 

0 

1 

0 


0 


y = [l 9 20]x 

which only indicates a permutation in the state space variables, 
that is 


X = 


0 

0 

1 


0 

1 

0 


1 

0 

0 


X 


Employing the partial fraction expansion (which can be obtained 
by the MATLAB function residue), the transfer function is 
written as 

Y(s) (5 + 5)(5 + 4) 6 6 1 

U(s) (5 + 1)(5 + 2)(5 + 3) 5 + 1 5 + 2~^5 + 3 
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The state space model, directly written using (3.29), is 



"-1 

0 

0 ■ 


"l" 

X = 

0 

-2 

0 

x + 

1 


0 

0 

-3 


1 


y = 

[6 

-6 

l]x 



Note that the parallel programming technique presented is valid 
only for the case of real distinct roots. If complex conjugate 
roots appear they should be combined in pairs corresponding to 
the second-order transfer functions, which can be independently 
implemented as demonstrated in the next example. 

Example 3.5: Let a transfer function containing a pair of 
complex conjugate roots be given by 


G{s) 


4 


-T + 


4 


5 + 1— J ' 5 + 1+J ' 5 + 5~^5 + 10 


: + 


We first group the complex conjugate poles in a second-order 
transfer function, that is 


, . 85 + 8 2 3 

G(s) = - 7 , -^-^- 

Then, distinct real poles are implemented like in the case of parallel 
programming. A second-order transfer function, corresponding 
to the pair of complex conjugate poles, is implemented using 
direct programming, and added in parallel to the first-order transfer 
functions corresponding to the real poles. The simulation diagram 
is given in Figure 3.5, where the controller canonical form is used 
to represent a second-order transfer function corresponding to the 
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complex conjugate poles. From this simulation diagram we have 

xi = —5xi + u 

X2 = —10X2 + U 
X^ = Xi 

Xi = — 20^3 — 20^4 + u 
y = 2 xi + 30^2 + 80:3 + 80:4 

so that the required state space form is 


"-5 

0 

0 

0 ■ 


"l" 

0 

-10 

0 

0 

x + 

1 

0 

0 

0 

1 

0 

0 

0 

-2 

-2 


1 


y=[2 3 8 8 ]x 



Figure 3.5: Simulation diagram for a 
system with complex conjugate poles 


o 




114 


STATE SPACE APPROACH 


Multiple Real Roots 

When the transfer function has multiple real poles, it is not 
possible to represent the system in uncoupled form. Assume that 
a real pole pi of the transfer function has multiplicity r and that 
the other poles are real and distinct, that is 

y(5) _ N{s) 

U{s) {s+piJ{s+PrJ,l)'"{s+Pn) 

The partial fraction form of the above expression is 

y(^) _ fell ^ ki2 ^ ^ hr ^ fer+1 ^ ^ fen 

U{s) S-\-pi (^s+pif {s+piY S-\-pr+l S-\-pn 

The simulation diagram for such a system is shown in Figure 3.6. 



Figure 3.6: The simulation diagram for the Jordan canonical form 
Taking for the state variables the outputs of integrators, the state 
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Space 

model is obtained 

as follows 






-pi 1 

0 

... ... 

0 

0 

... ... 

0 


0 -pi 

1 

0 

0 

0 

... ... 

0 


0 

0 

-pi 1 

0 

: 

: : 

: 

A = 

0 

. . . 

0 -pi 

1 

0 

... ... 

0 


0 

. . . 

0 

-Pl 

0 

... ... 

0 


0 0 

. . . 

. 

0 

'Pr+1 

0 

0 


0 0 

. . . 

. 

0 

• , 

—Pr+2 

i 


• • 

• 

• • 

• 

• 


0 


0 0 

. . . 

. 

0 

0 

0 

-Pn 







(3.26) 



B = [0 0 

. . . 

••• 0 1 

1 ••• 

. . . 

if 


—1 * 

. . 

• • ki2 hi 

h+1 

kr+2 

• • • kn], 

D = 


This form of the system model is known as the Jordan canonical 
form. The complete analysis of the Jordan canonical form requires 
a lot of space and time. However, understanding the Jordan form 
is crucial for correct interpretation of system stability, hence in the 
following chapter, the Jordan form will be completely explained. 

Example 3.6: Find the state space model from the transfer 
function using the Jordan canonical form 

5 ^ + 65 + 8 

(5 + 1) {s 3) 

This transfer function can be expanded as 

, 1.25 1.5 0.25 

G(s) =-^-^- 

s -\- 1 (5 + 1) 5 + 3 
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SO that the required state space model is 


"-1 

1 

0 ■ 


"o' 

0 

-1 

0 

x + 

1 

0 

0 

-3 


1 


= 1.25 -0.25]x 


3.4 The System Characteristic Equation and 
Eigenvalues 

The characteristic equation is very important in the study of both 
linear time invariant continuous and discrete systems. No mat¬ 
ter what model type is considered (ordinary nth-order differential 
equation, state space or transfer function), the characteristic equa¬ 
tion always has the same form. 

If we start with a differential nth-order system model in the 
operator form 

+ cin-iP^ ^ + • • • + Clip + ao^y(t) 

= {bmp'^ + bm-ip'^~^ H- \-bip-\- bo)u(t) 

where the operator p is defined as 

(P 

=-r^ , z = 1, 2,..., n — 1 
dP 

and m < n, then the characteristic equation, according to the 
mathematical theory of linear differential equations (Boyce and 
DiPrima, 1992), is defined by 

s Oji—xs a\s no ~ 0 

Note that the operator p is replaced by the complex variable s 
playing the role of a derivative in the Laplace transform context. 
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In the state space variable approach we have seen from (3.54) 
that 

G(s) = C(sl - A)“^B + D = pj——rC[adj(sI - A)]B + D 

= ^ ^ {C[adj(sl - A)]B + |sl - A|D} 

The characteristic equation here is defined by 


si — A = 0 


(3.28) 


A third form of the characteristic equation is obtained in the 
context of the transfer function approach. The transfer function of 
a single-input single-output system is 


+ a-n—lS^ ^ + • • • + + <20 


(3.29) 


The characteristic equation in this case is obtained by equating 
the denominator of this expression to zero. Note that for multi- 
variable systems, the characteristic polynomial (obtained from the 
corresponding characteristic equation) appears in denominators of 
all entries of the matrix transfer function. 

No matter what form of the system model is considered, the 
characteristic equation is always the same. This is obvious from 
(3.95) and (3.97), but is not so clear from (3.96). It is left as an 
exercise to the reader to show that (3.95) and (3.96) are identical 
(Problem 3.30). 

The eigenvalues are defined in linear algebra as scalars. A, 
satisfying (Fraleigh and Beauregard, 1990) 


Av = Av 


(3.30) 
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where the vectors v 7 ^ 0 are called the eigenvectors. This system 
of n linear algebraic equations (A is fixed) has a solution v / 0 
if and only if 

|AI-A|=0 (3.31) 

Obviously, (3.96) and (3.99) have the same form. Since (3.96) = 
(3.95), it follows that the last equation is the characteristic equa¬ 
tion, and hence the eigenvalues are the zeros of the characteristic 
equation. For the characteristic equation of order n, the number 
of eigenvalues is equal to n. Thus, the roots of the characteristic 
equation in the state space context are the eigenvalues of the ma¬ 
trix A. These roots in the transfer function context are called the 
system poles, according to the mathematical tools for analysis of 
these systems—the complex variable methods. 

Similarity Transformation 

We have pointed out before that a system modeled by the state 
space technique may have many state space forms. Here, we 
establish a relationship among those state space forms by using 
a linear transformation known as the similarity transformation. 

For a given system 

X = Ax + Bu, x(0) = Xo 
y = Cx + Du 

we can introduce a new state vector x by a linear coordinate 
transformation as follows 

x = Px 

where P is some nonsingular n x n matrix. A new state space 
model is obtained as 

X = Ax + Bu, x(0) = Xo 

y = Cx + Du 


(3.32) 
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where 


A = P ‘ AP, B = P 'B, C = CP, D = D, x(0) = P 4(0) 

(3.33) 

This transformation is known in the literature as the similarity 
transformation. It plays an important role in linear control system 
theory and practice. 

Very important features of this transformation are that under 
similarity transformation both the system eigenvalues and the sys¬ 
tem transfer function are invariant. 

Eigenvalue Invariance 

A new state space model obtained by the similarity transforma¬ 
tion does not change internal structure of the model, that is, the 
eigenvalues of the system remain the same. This can be shown 
as follows 


si-A 


sI-p-^AP 


P“^(sl-A)P 

Isl-Al 


• -1 


5l-A P 


(3.34) 


Note that in this proof the following properties of the matrix 
determinant have been used 


dot (Ml M 2 M 3 ) = detMi xdetM 2 XdetM 3 


see Appendix C. 


cletM^' 


1 

detM 


Transfer Function Invariance 

Another important feature of the similarity transformation is that 
the transfer function remains the same for both models, which can 
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be shown as follows 


G(s) = C - a) ^B + D = CP(sI-P^'AP) ^P^‘B + D 

= CP[P \sl-A)P]“^P 'B + D 
= CPP“^(sI - A) 'PP“^B + D 
= C(sl - A)^'B + D = G(s) 

(3.35) 

Note that we have used in (3.103) the matrix inversion property 
(Appendix C) 


(MiM 2M3)^‘ = Mg 

The above result is quite logical—the system preserves its in¬ 
put-output behavior no matter how it is mathematically described. 

Modal Transformation 

One of the most interesting similarity transformations is the one 
that puts matrix A into diagonal form. Assume that P = V = 
[vi, V 2 , ..., v„], where v* are the eigenvectors. We then have 

V~^AV = A = A = diag(Xi, A 2 ,..., A^) (3.36) 

It is easy to show that the elements A*, z = 1, on the 
matrix diagonal of A are the roots of the characteristic equation 
I^I — A| = I^I — A| =0, i.e. they are the eigenvalues. This can 
be shown in a straightforward way 

I^I — A| = det{diag(s — Xi,s — A 2 ,..., 5 — A^)} 

= (5 - Ai )(5 - A 2 ) • • • (5 - A„) 

The state transformation (3.104) is known as the modal transfor¬ 
mation. 
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Note that the pure diagonal state space form defined in (3.104) 
can be obtained only in the following three cases. 

1. The system matrix has distinct eigenvalues, namely Ai / A 2 / 

• • • 7^ A^. 

2. The system matrix is symmetric (see Appendix C). 

3. The system minimal polynomial does not contain multiple 
eigenvalues. For the definition of the minimal polynomial and 
the corresponding pure diagonal Jordan form, see Subsection 

4.2.4. 

In the above three cases we say that the system matrix is diago- 
nalizable. 

Remark: Relation (3.104) may be represented in another form, 
that is 

or 

W^A = AW^ 

where 

= V“‘ ^ W^V = I 

In this case the left eigenvectors w*, z = 1, 2,n, can be computed 
from 

rji 'T' 'T' 

A = A*w^ A w* = A*w* 

where W = [wi, W 2 ,..., w„]. Since |AI — A| = AI — AJ , then 
A* is also an eigenvalue of AJ . 

There are numerous program packages available to compute 
both the eigenvalues and eigenvectors of a matrix. In MATLAB 
this is done by using the function eig. 
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3.4.1 Multiple Eigenvalues 

If the matrix A has multiple eigenvalues, it is possible to transform 
it into a block diagonal form by using the transformation 

J = V^^AV (3.37) 


where the matrix V is composed of n linearly independent, so- 
called generalized eigenvectors and J is known as the Jordan 
canonical form. This block diagonal form contains simple Jordan 
blocks on the diagonal. Simple Jordan blocks have the given 
eigenvalue on the main diagonal, ones above the main diagonal 
with all other elements equal to zero. For example, a simple Jordan 
block of order four is given by 

A^ 1 0 0 ■ 

0 A^ 1 0 

0 0 A^ 1 

0 0 0 A^ 

Let the eigenvalue Ai have multiplicity of order 3 in addition 
to two real and distinct eigenvalues, A 2 / A 3 ; then a J matrix of 
order 5x5 may contain the following three simple Jordan blocks 


Jz(A^) = 


J = 


Ai 

0 

0 

0 

0 


10 0 0 
Ai 1 0 0 

0 Ai 0 0 

0 0 A 2 0 

0 0 0 A3 


However, other choices are also possible. For example, we may 
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have the following distribution of simple Jordan blocks 

"Ai 1 0 0 0 1 [Ai 0 0 0 0 " 

OAiOOO OAiOOO 

J=0 0Ai0 0 orJ=0 0Ai0 0 

000 A2O 000 A2O 

0000 A3 J [0000 A3 

The study of the Jordan form is quite complex. Much more about 
the Jordan form will be presented in Chapter 4, where we study 
system stability. 

3.4.2 Modal Decomposition 

Diagonalization of matrix A using transformation x = Vx makes 
the system x = Ax + Bu diagonal, that is 

= Ax + (V“^B)u = Ax + (w^b)u, x(0) = Xo 

In such a case the homogeneous equation x = Ax, x(0) = Xq, 
becomes 

X = Ax, x(0) = V“^x(0) = V“^Xo 
or 

OCi — A^^/^, Z — f, * * * , Tlj 

This system is represented by n independent differential equations. 
The modal response to the initial condition is 

±{t) = = e'^'W^Xo 

\t 


or 


x^{t) = = (^wfxo^e 
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The response x(t) is a combination of the modal components 

x(t) = Vx(t) = Ve^^V^o = Ve^^W^Xo 
= (^wfxo^e^^Vi + (w2^o^e^^^V2 H-h 

(3.38) 

This equation represents the modal decomposition of x(t) and it 
shows that the total response consists of a sum of responses of all 
individual modes. Note that w^^Xq are scalars. 

It is customary to call the reciprocals of A^ the system time 
constants and denote them by T^, that is 

'T~t - T 5 Tj - 1 5 2 5 ,^ Tl 

Az 

This has physical meaning since the system dynamics is determined 
by its time constants and these do appear in the system response 
in the form 

The transient response of the system may be influenced differ¬ 
ently by different modes, depending of the eigenvalues Some 
modes may decay faster than the others. Some modes might be 
dominant in the system response. These cases will be illustrated 
in Chapter 6. 

Remark: A similarity transformation A = AV can be used 
for the state transition matrix calculation. Recall 

x(t) = e^^x(O), x(t) = Vx(t), x(0) = Vx(0) 

and 

x(t) = V“^e''‘Vx(0) = 'j>(i)x(0) 

Hence 

$(t) = = W^e'^'V 
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or, in the complex domain 


<J>(s) = V“^(sl-A)^"V 


- 1 - 


V ^diag{s — Xi, s — X 2 ,..., s — Xn} 


-1- 


= V ^diag< 


S-Xl^ S-X2""' S-Xn 


>V 


Remark: The presented theory about the system characteristic 
equation, eigenvalues, eigenvectors, similarity and modal trans¬ 
formations can be applied directly to discrete-time linear systems 
with Ad replacing A. 



